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Abstract 

The role of the relativistic ampHtudes for a number of 0{k) processes usually 
neglected in potential model calculations of NN-bremsstrahlung is investi- 
gated. In particular, we consider the A-excitation pole contributions related 
to the one-pion and one-rho exchange and in addition include the exchange 
contributions induced by the radiative to, p ^ ttj decays. The contributions 
are calculated from relativistic Born amplitudes fitted to A-production and 
absorption data in the energy range up to 1 GeV and then used to sup- 
plement potential model and soft photon calculations for nucleon-nucleon 
bremsstrahlung. The effects on A^A^7-observables, although moderate in gen- 
eral, are found to be important in some kinematic domains. 
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I. INTRODUCTION 



Nucleon-nucleon bremsstrahlung has been extensively investigated both experimentally 
and theoretically during the past 30 years. For pp — »• pp'j the experimental data available 
covers an energy domain between 42 MeV and 730 MeV 0], the most recent results 
being obtained at incident proton energies around the pion production threshold Up 
to the pion production threshold, the NN'j potential model using realistic NN-potentials 
for the description of the nuclear force plus first order approximation for the electromagnetic 
interaction gives the most successful description. 

A comparison of different NN-interactions yields widely equivalent A^A^7-results in the 
whole kinematic range [§, so that bremsstrahlung provides a sensitive test for the dynam- 
ical model used to describe the photon emission. In the framework of the potential model, 
however, there is little controversy about the basic features and the correction terms to be 
used. In particular, relativistic spin corrections as derived in [0,^ and rescattering contri- 
butions [^] are used in the most recent analyses [§|,|10|) see also |Tl|jr^ . Moreover, the role of 



Coulomb corrections in pp7 [|13| and exchange currents in the np'j cross section have been 
discussed P,P . 



The main theoretical problems left are presumably related to shortcomings in principle 
of the potential model description. For example, Lorentz invariance can only be accounted 
for approximately by a covariant treatment of the kinematic transformations, inclusion of 
higher order terms in (p/m) in the electromagnetic operator (relativistic spin corrections) 
and an appropriate description of the NN-interaction. Further common features of existing 
potential model calculations are the absence of dynamical baryon resonances and the neglect 
of two-body currents beyond the 0{k^) terms given by the soft photon approximation (SPA). 

The role of the A-excitation was studied in the 1970's in two different approaches. The 
authors of used a dispersion analysis in order to correct the one-pion exchange neutron- 
proton electromagnetic current. The effect of the A-resonance on the nprf coss section at 200 
MeV turned out to be small. In |T6[ and |]T^, the A-excitation part of the pp^ amplitude was 
derived from phenomenological Lagrangians and combined with soft photon or one-boson 
exchange Born calculations for the radiative background. 

A more sophisticated potential approach including effects of the A based on a coupled 
channel calculation of the half off-shell NN- and NA T-matrices together with a phenomeno- 



logical A^A7- vertex has been published recently |T8,19|. In contrast to the ?T,p7-results of 



15| , the authors find an appreciable A contribution to the 280 MeV pp^y observables. 

For higher energies where the potential model is inappropriate, the closely related prob- 
lem of dilepton production has recently been studied in the framework of an effective one- 
boson exchange plus A-excitation Born approximation [pp]]. 

As mentioned above, a rigorous derivation of induced two-body currents is not possible 
to all orders. The reason is that, unlike e. g. for a one-pion exchange potential the gauge 
invariant replacement Vn{p) ^ Vm (p — eAj in the argument of a general NN-potential 
PJT^ leads to an unique expression for the induced current only in SPA. To this order, 



the induced current approximates contributions due to the exchange of charged mesons 
and vanishes for pp'j. The radiative decay processes we are considering here are thus not 
accounted for in the conventional potential model. So far, only the radiative uj — > tt^j decay 
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for J9J97 has been considered as an example of such an 0{k) internal radiation process p2[ . 



The purpose of the present work is to provide an estimate of the role of baryon reso- 
nances and internal radiation processes reliable up to photon energies of about 300 MeV 
corresponding to the highest photon energy possible in the 730 MeV j9p7-experiment of 
0. From a comparison with recent analyses of pion-photoproduction ||2^, we expect the 



A(1232)-resonance and the radiative decay of the u and p to be the leading corrections. 

A coupled channel calculation for the A^A^ — A^A system omitting the contribution of 
the AA-states p4[ is typically represented by a set of Lippmann-Schwinger equations 



Tnn — Vnn + VnnGnTnn + VnaGaTan, (1-1) 
Tan = Van + VanGnTnn + VaaGaTan- (1-2) 

Inclusion of the A thus modifies the NN-amplitude T/vat and yields an additional amplitude 
Tan- The pure NN-interaction below the pion production threshold is known to be well 
described by phenomenological and/or meson theoretic potentials so that we feel safe iden- 
tifying Tnn of equation ( |1 . 1|) with a pure NN-T- matrix calculated from a realistic nuclear 
potential in effect absorbing the VnaGaTan piece. By then adding the A-excitation Born 
terms, we basically neglect the modifications of the amplitude Tan (equation ( |1.2|) ) induced 
by the iterative terms. Empirically, i. e. from a comparison of the Born amplitude Van with 
coupled channel A-absorption predictions for lab energies of 50 — 300 MeV p^, and exper- 



imental A-production data at 800 MeV |25| and 970 MeV [^], it turns out that the effect 
of the iteration can be simulated by a simple energy-dependent rescaling of the Born ampli- 
tude in good approximation. We are then left basically with possible errors originating from 
the superposition of a (necessarily real and therefore not unitary) Born amplitude with the 
iterated potential model and SPA amplitudes. Our approach for the A is thus complemen- 
tary to a full coupled channel calculation []TB[ which systematically avoids the shortcomings 
mentioned above but cannot be extended easily to processes such as the internal radiative 
decays and - at least in the NNj-case - cannot be applied at higher energies. 

In the following two sections, a description of the model together with a discussion of 
the parameters will be given. Results for pp'j and np'-f are presented in section IV. We have 
taken care to test thoroughly the influence of the experimental and theoretical uncertainties 
for the various ingredients of the amplitude that can not be fixed to accurate data. The 
results might be useful to estimate the reliability not only of this but of any comparable 
calculation. 

II. INCLUSION OF THE A 

For the derivation of the A-excitation terms depicted in Fig. 1 we start from the inter- 



action Lagrangians for the isovector mesons vr and p (see ||2^) in the notation of |28 





WttNN 


^pNN — 


~9pNN 




9-kNA 
■) 







(2.1) 
(2.2) 

(2.3) 
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Cp^^ = t^^P^'^TrAid^Pu - d,p^) + H.c. , (2.4) 

TTlp 

— * 

where ipA is the spin-3/2 Rarita-Schwinger field and T stands for the isospin operator for 
the transition of an isospin-1/2 and isospin-1 to an isospin-3/2 particle |27]. We have used 
the pseudscalar vrA^A^-coupling for convenience. As the ttA^A^- vert ices in our calculations 
involve only on-shell nucleons, pseudovector coupling for the pion instead of equation 
would not change any of our results. From eq. - (p.4|), the vertex functions read 



K.M = K'.Ai) = iAg'" - 7V) 75 (2.5) 

^ nip 



and one derives the transition amplitudes for diagram (a) in Fig. 1 ||T] 

T('Hpi,P2;P3,P4) = u{p,)K,MPuAPi - k)epT';i^{k,p^ - k)uip,) 

+u{p,)K';iMPt{Pi - k)e^m^{k,pi - k)u{p,) 

Pl{q)u{p,Wp,,u{p,), (2.6) 

and equivalently for diagram (b) 

T('nPi,P2;P3,P4) = u{p;)e,v';i^{k,p^ + k)P^^{p^ + k)K,M<Pi) 

P"" {q)u{p4)A^NNu{p2) 
+uiPs)ep^'',,ik,p, + k)P^^{p, + k)K%^{q)u{pr) 

P^MHpM;.MP2). (2.7) 

Here, the energies and momenta are constrained according to pi+p2 = Ps+P^ + k, q = P4—P2 
is the four momentum transferred by the meson, the photon unit polarization vector, u, u 
denote free nucleon Dirac spinors and P the appropriate propagators for the mesons and 
the A. With the tilde in equation ( |2.7| ) we indicate that the vertex function has to be taken 
from the Hermitian conjugate of the corresponding Lagrangian. 

Associated with eqs. ( |2.6| ) and ( |2.7| ) are isospin factors, e. g. in case of proton-proton 
bremsstrahlung: 

2 

{Xp I ^3 I Xa+){Xa+ I ^3 I Xp){Xp \T3\Xp) = ^ 

The full A excitation amplitude including the isospin factors and a relative minus for graphs 
involving interchange of the final state particles is then 

Tr = |E{T«(P1,P2;P3,P4) + T»(P2,Pi;P4,P3) 

i=l 

- {T'-'\puP2]P4,P3) +tW(p2,pi;p3,P4))} 

= |E{(-T«(Pi,P2;P3,P4) + T«(p2,Pi;P4,P3)) 

i=l 

+ i-lY'^ {~T^'\puP2;P4,P3) + r«(P2,Pi;P3,P4))} . (2.8) 
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For the electromagnetic A^A current of diagram (a) we use the parametrization 



29 




mi 



75 ■ 



(2.9) 



T^ff^{k,p) is equal to r^jvA(^,p) up to a relative minus sign in the term proportional to Gi. 
The couplings Gi and G2 are experimentally determined by a fit to the Ml+/i?l+ multipole 
pion-photoproduction cross section. Values for Gi and G2 range from Gi = 2.0, G2 = 0, 
the prediction of the vector dominance model, to Gi = 2.68, G2 = —1.84 |18,29|, recently 
quoted values being Gi 
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2.208, G2 = -0.556 _ 
The choice of the A-propagator requires some caution. In |T^, the Rarita-Schwinger 
form 



p2 — + iTrriA 



g 



1 ^ , 2p>- 
3 3 mi 



3m\ 



(2.10) 



is used where the replacement m/^ —>■ mi^ — iT/2 is made in order to account for inelasticities 
due to pion production. As we wish to use Pa in the far off-resonance region (p^ << mA), we 
have to take the energy dependence of the width F ^ r(g) into account. For our purpose, 
we rely on the parametrization of ||30| requiring F to vanish below the pion production 
threshold: 



r(g^ 



N 



0, g < 0, q^N 



r(g,r7v) = 2-i{q^NR/m^f / (1 + {q^NR/m^f), ql^ > 0. 



(2.11) 



Here, q.,rN denotes the maximum momentum in the vrA^-subsystem of the process A^A^ 
NNtt: 



ql^ = {snN - {rri^ - mNf){snN - {rri^ + mNf)/4:SnN, 



{y/s - m^] 



(2.12) 



a/s is the invariant energy in the A^A^-system and R, 7 are adjustable parameters. Values 
of 7 = O.nMeV, R = 0.81 lead to a resonance width F = 120 MeV at s^n = 1236 MeV 
30| . This value might be slightly modified by a variation of 7. It has been shown |^ that 
this parametrization is in good agreement with the experimental ^33 phase as well as with a 



A self energy calculation P4 |. 

A^A^ — A^A potential models differ in the choice of the A-coupling constants as well as in 
the definition of the meson propagators and vertex form factors. We use these models as a 
starting point for our calculation defining the A-excitation amplitude up to a normalization. 
In a reliable calculation, the final results should not depend too much on the underlying pa- 
rameter sets provided they agree comparably well with the empirical NN-data. We therefore 



compare two different recent coupled channel models: The OBE-model of [2^|, denoted as 
model A, and a coupled channel version of the Bonn potential (model I of [^, p. 353), 
denoted as model B. For model A, we use the meson propagators 



mt 



pip) 

fJ,U 



-g^lv + q^qu/m 



mt 



(2.13) 
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The second term in P^^^ drops out in the calculations of this and the following section. 
The propagators of model B take the mass difference between the nucleon and the A into 
account: 



(2.14) 



In writing down the Lorentz invariant propagators ( p. 13 ) and ( 2.14 ) we have dropped the 
static approximation — of the potential models A and B. As we are far from the 

pole in all our geometries and rescale our final amplitudes this has very little effect on the 
end results. 

For the calculations using model A, each NNn/NNp-veitex is multiplied by a monopole 



cutoff Al/p/ 



/{A^/p - q ) and a dipole cutoff A^^/^/(A^^/^ + q^) applies at the NAn/p- 
vertices. Here again, we use Lorentz invariant instead of static expressions. For model B, 
both the NN- and A^A-vertices are regularized by cutoffs (A^ — rn^ip)/ (A^ — q^) where 
n = 1, 2 for TT, p exchange respectively. 

In general, the relativistic off-shell A-particle is allowed to propagate with both spin 3/2 
and spin 1/2. This extra freedom is reflected by additional off-shell terms in the interaction 
Lagrangians; the more general chiral vrA^A-vertex functions corresponding to the propagator 
( CT ) read |2|,g: 



A^ 



giTNA 



pNA 



(2.15) 



A generalization similar to A^^va (with parameter Z^) applies for the (leading) Gi-term 
of equation ( |2.9| ). The Z-parameters are not well determined either theoretically or ex- 
perimentally. Whereas the simple coupling scheme corresponding to Z = —1/2 is gener- 
ally preferred in potential models, Olssen and Osypowski find experimental values of 
Zj^ = 0±l/4, = l/4±l/4. Note that the off-shell parameters do not affect a cou- 
pled channel calculation at the OBE level but directly enter the Born amplitude for the 
bremsstrahlung process where the A is off-shell. We therefore rely on the parameters of 
Table 1 as fixed by NN-data (to all orders of a coupled channel calculation) but check the 
sensitivity of our results to a variation of Z. 

This completes the necessary ingredients for the A-excitation part of the bremsstrah- 
lung amplitude ( p.8| ). With the definitions of this chapter we have tried to incorporate as 
much as possible of the experimental information on the A^A excitation channels in the 
NN-interaction which have become available since the earlier Born calculations [p!6[ . 
On the other hand, the relativistic approach allows us - in contrast to a pure potential 
model calculation 



TBI - to extend the analysis to energies far beyond the pion production 



threshold. 
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TABLES 



Model A [121 



Model B [571 



4-JT 



A[MeV] 



9_ 

4lT 



A[MeV] 



NNtt 
NNp 

NAtt 
NAp 



14.16 
0.43 
K = 5.1 
0.35 
4.0 



1140 
1140 

910 
910 



14.4 
0.7 
K = 6.1 
0.35 
19.0 



1800 
2200 

920 
1140 



TABLE 
etc. see text 



L Meson parameter of model A and B. For the definition of form factors, propagators 
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III. INCLUSION OF THE RADIATIVE w/p DECAY 

We write the Lagrangian for the couphng of the cu-meson to the nucleons in analogy to 
equation (|2.2| ): 

i^uNN = -9u (^l^ipuj^ + ■£^'4!a^"'ij{d^uj^ - d^ujf,)^ , (3.1) 

and parametrize the radiative decay vertices as in p3| 

jC^TT-y = gu.n'yeupasid'' AP){d''iT°)uj^ -> A^^pPi^, q) = ~igu>-K^eyp„sk''q'' , (3.2) 

C-p^^y = gp.^^e.p„s{d''A<l^Wn')p'^'- ^ = 0, +, -. (3.3) 

With the meson propagators as in ( |2.13D and a vertex function A^^^^ as in eq. ( |2.5| ) this 
yields for the diagram (c) of Fig. 1 

Ti;VPl'P2;P3,P4) = n(P3)A;:^^(p3 -Pl)«(Pl)<)(P3 -Pl)A5;"(fc,P4 -P2)e^; 

X ^w(P4 -^2)^(^4)75^(^2) , (3.4) 
and consequently (see equation ( |2.8|) ) 

- {T{P1,P2;P4,P3) +T{p2,Pi;P3,P4)), 

T^^ = T{pi,p2;p3,P4) -T{p2,pi;p4,P3). (3.5) 

The p-decay amplitudes are obtained from T^^^ by interchange of the masses and cou- 
pling constants and multiplication with the isospin operator ti ■ T2. The matrix elements of 
this operator yields factors of 1 and for the neutral and charged pion decay amplitude in 
ppy. For npy the corresponding values are -1 and 2. 

The coupling constants g^T^^ and (^^7^7 are determined from the experimental decay widths 
of the vector mesons ||34|| : 

( 2 _™2\3 f 716 ± 75 keV for cj, 
T{uj/p ^ 7r7) = ^"/""^^ ^ = 121 ± 31 keV for p", (3.6) 



'^/p y 68 ± 7 keV for p±. 



This leads to the numerical values 



222 

^ = 0.715 ■ 10-=^m;2; ^ 0.125 ■ lQ-^m~^- ^ = 0.070 ■ lO-^m^l (3.7) 

An An An n \ j 

For the remaining coupling constants and cutoffs, we again follow the rationale of 
the previous section and use a complete set that has been successfully tested in pion- 
photoproduction ||23|. Thus we put 



^ = 0.563; Kp = 3.71; ^ = 5.07; k<^ = -0.12. (3.8) 
An An 

These values are in agreement with the quark model and the vector dominance assump- 
tion. All form factors in this section are set to 1; the pion couplings are taken from table 
I. Note that a measurement of the vector meson radiative decay contributions in the pion 
photoproduction process determines the product of the couplings in eqs. (|3.7|) and (^.81). 
With the well-known empirical value of gj^, the amplitude (p. 41) is therefore essentially fixed. 
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IV. RESULTS 



As mentioned earlier, the amplitudes ( |2.8| ) and ( |3.4| ) must be divided by an energy 
dependent scaling factor g{Tiab) in order to compensate for neglecting higher order terms in 
the Born series. From a comparison of our Born results using model A with the total A- 
absorption cross sections of calculated in the iterated coupled channel model, one finds 
a factor g = 1.65 at lab energies below 300 MeV. A calculation of the total A-production 
cross section and comparison with the experimental data at 800 MeV |^ and 970 MeV 
26| suggests a weak energy dependence of g so that 5'(730 MeV) 

m 



1.4. The latter value 
As we have to compensate also for the the 
1|), i^A^) , we use g{280 



is in close agreement with the result of 
change from static to Lorentz invariant propagators in eqs. 
MeV) = 1.95 and 51 (730 MeV) = 1.65 respectively. 

If we suppose the higher order corrections to be of equal importance for both model A 
and B and for the a;/p-decay amplitudes for which no experimental comparison could be 
made, the model is completely determined with the empirical values of g{Tiab) found above. 
The rescaled amplitudes ( |2.8|) and (|3.4| ) can thus be coherently added to the corresponding 
transition matrices of the potential model and the soft-photon approximation p5| , |36| so 
that the effects upon both cross sections and spin observables can be studied. The relevant 
formula for the observables considered can be found in the literature ( [^|I7|). For the 
potential model calculations we shall focus on the domain around pion production threshold 
where both correction effects should be maximum; the soft photon model will be used to 
reanalyze the experimental data of . 

In principle, with all conventions carefully chosen, the formulas in the last two sections 
give the correct relative signs for the corrections with respect to the leading (potential model 
or SPA) amplitudes. For an independent check we combined our amplitudes ( p.8|) and ( p.4|) 
with the pure one-pion exchange bremsstrahlung amplitude (Fig. 1, diag. (a) with the in- 
ternal A-line replaced by a nucleon line). If the exchanged pion is numerically put on-shell 
(g° — > + ml), this is equal up to a common factor to the pion-photoproduction ampli- 
tude with and without intermediate A excitation. In a near-threshold geometry (g^ ~ 0), 
the A-amplitude gives enhancement of the cross section by a few percent, in agreement with 
Peccei's chirally invariant Lagrangian for n^p p~f whereas the a;/p-decay amplitudes 
interfere destructively, as in [l23| . 

The sign of the u;-coupling constant given in eq. ( p.7| ) is consistent with the pion- 
photoproduction data ||2^. The sign of the much smaller radiative p-decay amplitude, 
however, is still the object of controversies [^. We adopt here the sign convention of Gari 
et al. where the p-decay enhances the contribution of the cj-decay in ppy; as the p- 



contribution is only about 2 % of the ui (compare the respective coupling constants), a switch 
in the sign of gpT,y would not alter our conclusions. 



A. Potential Model Results 

All the calculations in this section are based on the inversion potential to the Nijmegen-II 
NN-phase shifts used in Ref. 0. Using another realistic NN-potential would not change 
any of our results. The basic features of the potential model can be essentially represented 
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in terms of two variables: the total energy of the process Ti^b cis a crucial parameter for the 
relative size of the various contributions in the amplitude, and the opening angles of the 
two protons determining the maximum photon energy and thus the off-shell signature of 
the process [^. The two examples of Fig. 2 show the coplanar, exclusive pp7 cross section 
d'^cr/dflidfl2d9^ and the analyzing power Ay for the smallest and largest outgoing proton 
angle pairs measured in the Tiab = 280 MeV TRIUMF experiment 0]. The plots show that 
both A-excitations and internal radiative decay contributions are small for the cross section 
at small proton angles as well as for the analyzing power Ay at large proton angles. The 
A-contribution is the leading effect but tends to be partly cancelled by the radiative decay 
contributions. The maximum net effect of both contributions can increase the differential 
cross section by 15 — 20% at medium photon angles but amounts to only a few percent in 
the total cross section. 

Among the examples of Fig. 2, only the 12.4°, 14° analyzing power shows a sizeable 
dependence on the parametrization of the vr and p-exchange. The choice of different elec- 
tromagnetic couplings, e. g. Gi = 2.68, G2 = —1.84, increases the A-contribution in the 
cross sections by about 20 % which translates to an enhancement of the complete ppj cross 
section by about 5 % at maximum (Figure 3). Conversely, the choice of Z^^ = = —1/4 
PB| reduces the size of the A-contributions by roughly 20%. These two choices for the 
couplings and off-shell parameters are likely to be rather extreme. Moreover, because of its 
small relative size, a rescaling error in the cu/p-amplitude would not affect much the end 
results. One might therefore conclude that the corrections are reasonably well determined 
in our model. 

In order to obtain a more general picture we calculate the double differential cross sec- 
tion da/dQidQ2 = J d9^{da / dVLidVL2d9^) and represent its relative enhancement due to the 
correction terms as a function of the symmetric angle 6 = 61 = 62 (Figure 4). Note that 
small 6 values correspond to a suppression of the meson four momentum transfer {q —>■ 0) 
whereas for large 6, the elastic limit ^ is reached. In both cases, the amplitudes ( ^.8| ) 
and ( |3.4|) are suppressed: they are maximum in the medium proton angle region 6 ~ 20°. 
The relative size of the corrections dies down as the energy decreases and reaches only about 
3% at Tiab = 200 MeV. 

In Figure 5 the analysis of Fig. 2 is repeated for the np7-observables. Here, the isospin 
factors yield strong destructive interference between the various Feynman graphs so that 
both the A excitation and the radiative decays become negligible corrections to the potential 
model cross section. The same result holds for the np7-analyzing power. We thus confirm 
the result of that A-corrections to the n]97-amplitude are weak. 



B. 730 MeV (SPA) Results 



Two typical geometries of the 730 MeV pp7-experiment reported in Ref. |3^ are shown 
in Fig 5. One proton is emitted at 61 = 50.5°, 0i = 0° and the polar angles of the photon 



are 9^ 



179° for counter G7 and 9^ = 54°, 



131° for counter GIO. The 



solid curves have been obtained with the soft photon approximation of [^. Contributions 
in 0{k) are only partly included in SPA; the corrections considered here are thus specific 
candidates for the missing 0{k) effects. 
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As expected the corrections are negligible for small k but increase with the photon 
energy. Analogous to the 280 MeV examples, there is cancellation between the A- and 
cij//7-effects. The actual size is again geometry dependent but is limited to about 20% of the 
SPA cross section at k = 150 MeV. A variation of the electromagnetic coupling constants 
and off-shell parameters within the experimental limits yields effects analogous to those 
shown in Fig. 3 and leaves a freedom of a few percent in our final k = 150 MeV results. 
Correspondingly, the results are not much altered by interchange of model A and B. We 
conclude that a A-amplitude consistent with experimental results can only resolve a small 
part of the discrepancy between the SPA and the 730 MeV pp7-results at photon energies 
above ~ 100 MeV. 



V. CONCLUDING REMARKS 

We have evaluated the A-excitation and radiative cu/p-decay corrections to proton- 
proton and neutron-proton bremsstrahlung for the energy range up to Tiab ~ 1 GeV. This 
was done by calculating the relativistic Born amplitudes and adding them to potential model 
and SPA amplitudes. The Born amplitudes were normalized by calculating A-production 
and absorption cross sections in the same model and fitting to experimental data and coupled 
channel predictions. 

Of the two processes considered, the A-excitation dominates but is generally partly 
compensated by the radiative decay contributions. Both corrections together increase the 
280 MeV ppj integrated cross section (Pa / dVLidVL2 by an amount depending essentially on 
the proton opening angles 9i and reaching a maximum of roughly 7.5% at 6i 62 ^ 20°. 
The relative effect on d^a/dQidQ2dO^ is maximum for photon emission around 6^ = 90" and 
is typically < 20%. 

In the 730 MeV Rochester geometry, the corrections become relevant for photon energies 
above ~ 100 MeV but are not big enough to complement the soft photon approximation 
and fit the experimental data . 

We have taken care to study and discuss the limits and uncertainties of our model. The 
results have been obtained on the basis of two different parametrizations of the A-excitation 
(model A and B of Table 1) but show little sensitivity to the underlying coupled-channel 
model. The fact that the A-excitation cross sections of Model A and B are nearly identical, 
despite the rather massive differences in the values of the p-couplings and the choice of 
form factors and propagators, is reassuring for our approach but corresponds to the result 
obtained earlier that rather different potential models, once they fit the NN-scattering data, 
yield very similar bremsstrahlung results ^ and is disappointing if one hopes to understand 
the underlying physical reaction mechanism from NN^ measurements. 

The experimental uncertainty in the radiative A-decay constants and the off-shell free- 
dom of the A suggest a theoretical error of about ±20% in the final correction amplitudes, 
i. e. up to ~ ±4% in the total (280 MeV) pp'j cross section. We should mention that the 
discrepancies between the ppy predictions of different NN-potential models, which most 
bremsstrahlung experiments in the past were designed to isolate, are of the same order of 
magnitude. Given the present experimental status of nucleon-nucleon bremsstrahlung and 
the theoretical ambiguities mentioned, it is not possible to use A^A^7-data for putting limits 
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on the parameters entering our correction amplitudes. 

In view of the similarities of our calculation with previous ansatzes, a comprehensive 
comparison of the results seems worthwhile. First note that we have eliminated the un- 
certainty with respect to the relative sign of the relativistic corrections stated by Kamal 
and Szyjewicz through a comparison with the pion-photoproduction process. Our uj / p- 
decay amplitude extends the amplitude of ||22| by including the w-tensor coupling terms and 
the p-decays and is, if we divide by g{Tiab), slightly smaller than given there. 

The discrepancy of our result is more serious for the 730 MeV A-excitation calculation 
where the previous authors lack a reliable model for the NN-channel pp7-amplitude. The 



authors of therefore essentially fit the A-excitation part to the experimental pp7 data. 
In ||16|, the suppression of the amplitude which we simulate by appropriate form factors is 
neglected. In both cases, the pure A-cross section becomes bigger than ours by more than 
a factor of 10 and would be in clear contradiction to the 800 MeV A-production data of, e. 
We also stress the importance of taking interference terms correctly into account. 



The dispersion theoretic approach used in ||T5[ to estimate the role of the A-resonance 
in a 0PE-rap7 calculation is rather different from our model. Correspondingly, the results 
agree only in relative size. 

Finally, a comparison with the coupled-channel A-results of []TB[ in the low energy region 
shows good agreement for the cross sections and even for the analyzing powers Ay. This is 
encouraging for the feasibility of our method as well as for the various extensions made. 

We wish to thank Dr. S. Scherer for valuable discussions. This work was supported in 
part by a grant from the Natural Sciences and Engineering Research Council of Canada. 
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FIGURES 



FIG. 1. A-excitation and internal radiation processes considered in this text. 

FIG. 2. Coplanar ppj exclusive cross sections and analyzing powers at Tiab = 280 MeV and for 
the smallest and largest proton angle pairs di, 62 of the TRIUMF experiment Q (the datapoints 
shown are rescaled with a factor of 0.67 as in this reference). The curves denote the pure potential 
model of Ref. |^ (solid), potential model plus pure uj/p-dec&y (long dashed), and the full model 
according to model A (dotted) and B (dashed-dotted). 

FIG. 3. Typical uncertainties of the A-contribution. The curves show the pure potential 
model of Ref. |^ (solid) and the full calculation according to model A with standard parameters 
Gi = 2.208, G2 = —0.278, Z,^ = Zj = —1/2 (long dashed). The upper bound of the shaded area in 
the cross section and lower bound in Ay corresponds to larger electromagnetic couplings Gi = 2.68, 
G2 = —1.84. The lower bound in the cross section and upper bound in Ay is calculated with larger 
Z-parameters Zj^ = Z^ = —1/4. 

FIG. 4. Relative change of the integrated, coplanar pp-y cross section (f a / dQ,idQ2i^^^^ 
corrections)/d^(T/(iS7i(i02(without corrections) at 280 MeV (dashed), 200 MeV (dashed-dotted) 
and 157 MeV (dotted) for parameter set A as a function of the symmetric proton angles 61 = 62- 

FIG. 5. Coplanar npj exclusive cross section and analyzing power at Tiab = 280 MeV. The 
curves denote the pure potential model of Ref. |6| (solid) , plus pure lo / /3~decay (long-dashed) , and 
the full model according to model A and B (dashed-dotted and dotted). 

FIG. 6. 730 MeV pp'y measurements and soft photon approximation in two typical geometries 
(see (1979) and text). The curves denote the pure SPA results (solid), plus pure lo /p-decay 
(long-dashed), and the full model according to model A (dotted). 
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